import numpy as np
from zhnbodeplot import plot_bode
index = np.linspace(-1, 3, 41)
omega = np.power(10, index)

phi = -np.arctan2(omega, 3.3)*180/np.pi
amp = 11 * (0.3 + 10/omega**2)
amp /= np.sqrt(omega**2 + 3.3**2)
# phi -= np.arctan2(omega, 1)*180/np.pi
# amp /= np.sqrt(omega**2 + 1)
# phi -= np.arctan2(omega, 1)*180/np.pi
# amp /= np.sqrt(omega**2 + 1)
plot_bode(index, amp, phi, False)